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Abstract 

The magnetic response of irreversible type-II superconductor slabs subjected to in-plane rotating mag- 
netic field is investigated by applying the circular, elliptic, extended-elliptic, and rectangular flux-line-cutting 
critical-state models. Specifically, the models have been applied to explain experiments on a PbBi rotating 
disk in a fixed magnetic field H a , parallel to the flat surfaces. Here, we have exploited the equivalency of 
the experimental situation with that of a fixed disk under the action of a parallel magnetic field, rotating 
in the opposite sense. The effect of both the magnitude H a of the applied magnetic field and its angle of 
rotation a a upon the magnetization of the superconductor sample is analyzed. When H a is smaller than 
the penetration field Hp, the magnetization components, parallel and perpendicular to H a , oscillate with 
increasing the rotation angle. On the other hand, if the magnitude of the applied field, H a , is larger than 
Hp, both magnetization components become constant functions of a s at large rotation angles. The evolution 
of the magnetic induction profiles inside the superconductor is also studied. 

PACS numbers: 74.25.Ha, 74.25. Op, 74.25.Sv, 74.25.Wx 

Keywords: flux cutting, flux transport, vortex pinning, critical state, hard superconductor 
I. Introduction els in the past few years. 



The discovery of the phenomenon known as quasisym- 
metrical collapse of magnetization [T|, which is ob- 
served in superconductors subjected to crossed mag- 
netic fields and well interpreted within the simple 
Bean's critical-state model [2\ |3j, has been a turning 
point in the understanding of the magnetic behavior 
of hard (irreversible type-II) superconductors. Until 
then, the generalized double critical-state model (GD- 
CSM) m [8], which is based on fundamen- 
tal physical concepts such as flux transport and flux- 
line-cutting [5J [ID], was successfully employed to ex- 
plain a variety of experiments where flux cutting oc- 
curs [III OU ECU [H ED EE]- An important feature of 
the GDCSM is the assumption that flux cutting and 
flux depinning do not affect each other. Besides, the 
GDCSM is inherently anisotropic because the thresh- 
olds for these two effects are given by two independent 
parameters, namely the critical current densities par- 
allel (J c y) and perpendicular (J c j_) to the local mag- 
netic induction B. However, since the GDCSM cannot 
reproduce the features of magnetic moment collapse 
|17l US] , whereas isotropic Bean's model does it, the 
main assumption of the GDCSM has been questioned, 
motivating the development of new critical-state mod- 



In Ref. |19) . the so-called elliptic flux- line-cutting 
critical-state model was proposed. This model intro- 
duces the anisotropy, induced by flux-line-cutting ef- 
fects, by using a procedure similar to that for struc- 
turally anisotropic superconductors |20[ [21] . i.e. the 
magnitude of the critical current density J c . being 
the only parameter used within the isotropic Bean's 
model, is substituted by a symmetrical tensor (J c )ik 
with principal values </ c || and J c j_, corresponding to 
the directions along and across the local magnetic in- 
duction B. In good agreement with the experiment 
on YBa2Cu 3 07_a samples QUE], the elliptic critical- 
state model predicts the quasisymmetrical suppression 
of the average magnetization < M z >, for param- 
agnetic and diamagnetic initial states, by sweeping a 
transverse field H y of magnitude much smaller than 
dc-bias magnetic field H z |19l 122] . When the mag- 
nitudes of the crossed fields H y and H z are compa- 
rable, the value of the magnetization < M z > after 
many cycles of the transverse field H y turns out to be 
positive for both diamagnetic and paramagnetic initial 
states if J c u > J c ±. To our knowledge, the observa- 
tion of such a paramagnetism of hard superconduc- 
tors was first reported in Refs. [231 124|. The elliptic 
model also describes the behavior of < M y > (H y ) 
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the effect of the magnitude H a of the applied mag- 
netic field upon the dependencies of the magnetization 
components, parallel and perpendicular to H a , on the 
rotation angle of the superconductor disk. The evolu- 
tion of magnetic induction profiles will also be studied 
to explain the magnetic response of the rotating hard- 
superconductor sample. 



II. 



Theoretical formalism 



Let us consider a superconducting slab of thickness d, 
which occupies the space < x < d and is subjected 
to a magnetic field H a parallel to its surfaces: 

H a d. s = H a [ysm(a s ) + z cos(a s )], (2.1) 



H, 



and < M z > (H y ) in crossed fields H y and H z [19, 25j, the experiment [33] (Sec. III.). Here, we shall analyze 
which was observed in the experiments on a VTi ribbon 
with nonmagnetic initial state [b"fll26| . Here, the good 
agreement with the experiment was achieved by using 
a relatively large anisotropy parameter J c \\/J c ±_ — 6. It 
should be noticed that the Bean's critical state model 
predicts neither the phenomenon of the paramagnetism 
of hard superconductors nor the behavior of the com- 
ponents of the average magnetization found in Refs. 
[HUM]. Furthermore, as it is shown in Refs. |19[ I27|. 
the elliptic critical-state model successfully describes 
the magnetic response of superconducting disks under- 
going oscillations in a magnetic field of fixed magni- 
tude for nonmagnetic, paramagnetic, and diamagnetic 
initial states [TT] . 

Despite the great success of the elliptic model |19| . 
it turns out that there exist phenomena, associated 
with flux cutting, which are not completely described 
within such a model. So, in a very recent work [28j, 
the elliptic critical-state model and other four theo- 
retical approaches for describing the critical state of 
type-II superconductors (GDCSM, extended GDCSM 
[2UI 150] . extended elliptic critical-state model [251 131j . 
and an elliptic critical-state model based on the vari- 
ational principle [32]) were tested. There, the angular 
dependencies of the critical current density J c and the 
electric field E (for J just above J c ) were measured, 
using an epitaxially grown YBCO thin film, and com- 
pared with the predictions of the five theories. The 
measurements of angular dependence of the critical- 
current density J c demonstrated a behavior rather sim- 
ilar to that assumed by the elliptic critical-state mod- 
els. Besides, the smooth angular dependence of the 
ratio of the transverse to the longitudinal components 
of the electric field E y /E z for J just above J c , pre- 
dicted by the three elliptic models, was verified in the 
experiment [25]. However, the original critical-state 



where a s is the angle of the applied magnetic field H a 
with respect to the z-axis. Hence, the magnetic in- 
duction B(x, t) inside the superconducting slab can be 
expressed as 

B = B(x,t)[ysin(a(x,t)) + zcos(a(x,t))], (2.2) 

where B and a are respectively the magnitude and the 
tilt angle of the magnetic induction. It is convenient 
to write the electric field E(x, t) and the electrical cur- 
rent density J(x,t) in terms of their components par- 
allel and perpendicular to the local magnetic induction 
B(x,t): 



E(x,t) = 
J(x,t) = 



E\\{x,t)a{x,t) 
Jy (x, t)a(x, t) - 



-E ± (x,t)3{x,t),{2.3) 
J x {x,t)$(x,t),(2A) 



model [19J leads to small values of the ratio E y /E z in 
comparison with the experimental data and the results 
obtained from the other two elliptic models. On the 
basis of this detailed comparison between experiment 
and the five theories, it was concluded in Ref. [35] that 
the experiment favors only one of the models, namely 
the extended elliptic critical-state model. 

The aim of the present work is to investigate the 
behavior of a hard superconductor in a parallel rotat- 
ing magnetic field (or equivalently, the response of a 
rotating superconductor in a fixed magnetic field) and 
to compare the predictions of four critical-state mod- 
els with experiment. Concretely, we shall consider the 
Bean's critical-state model [3J [3J, the original elliptic 
critical-state model [191 122] . the recently-proposed ex- 
tended elliptic model (S5J[3T], as well as the GDCSM 
|H 13 El 13 H], whose main characteristics and assump- 
tions will be revisited in Sec. |IL] We shall numeri- 
cally solve Maxwell equations with the material equa- 
tion postulated by each of the considered critical-state 
models to calculate magnetization curves for a super- 
conductor disk rotating in a fixed magnetic field as in 



where (3(x,t) = x x 6c(x, t). Inside the supercon- 
ductor sample, we shall assume that the magnetic 
induction and the magnetic field satisfy the relation 
B(x,t) = fiQH(x,t), which is good enough for applied 
magnetic fields much larger than the first critical field 
(H a ^> H c i). Moreover, any surface barrier against 
the flux entry (or exit) will be neglected. According 
to the planar geometry of the problem, we can rewrite 
Ampere and Lorentz laws, 



as follow 



V x B(x,t) 

V x E(x,t) 

dB 

dx 
da 



B 



dx 



8Ej_ 
dx 



E 



E± 



da 
dx 



da 

l dx~ 
dE l 

dx 



dB 

-fM)J\\, 

dB 

~~dt' 
da 



-B 



dt 



(2.5) 
(2.6) 

(2.7) 
(2.8) 

(2.9) 
(2.10) 



To solve the resulting system of differential equations 
for E, B and J, one should add the material equa- 
tion. Below, we shall use the material equations corre- 
sponding to the circular, elliptic, extended-elliptic, and 
rectangular flux-line-cutting critical-state models. 
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A. Circular model 

The first model for describing the magnetic behavior 
of superconductors in multicomponent situations was 
proposed by Bean (2j [3] - According to it, the critical 
current density J points always along the local electric 
field E. Hence, 



J = Jr. 



E 



(2.11) 



The magnitude of the critical current density J = J c is 
the unique phenomenological parameter used and may 
depend on the magnitude of the magnetic induction 
B. In the planar geometry [see Eqs. (2.1 1-( 2.4 )] , the 



assumption J = J c corresponds to a circle in the J±- 
Jy plane. 

In numerically solving the system of equations 
( 2.7 )-( 2.10) for the electromagnetic fields, it is neces- 



sary to rewrite Eq. (2.11 1 as 



E(J) 



(2.12) 



E{J) 



(2.13) 



0, J < J C {B) 

p(J-J c (B)), J>J C (B) 

where p is an effective resistivity. It should be men- 
tioned that for slow variations of the surface bound- 
ary conditions, producing a small magnitude of the in- 
duced electric field (E <C pJ c ), the magnetic induction 
profiles are practically relaxed and independent of the 
parameter p |34| . 

B. Elliptic model 

The elliptic flux-line cutting critical-state model |19[ 
122] 125] postulates: 



Ji = (Jc) 



i k 



Ek 

E 



where 



{Jc)ik — Jc,i{B) Si 



,k=±,\ 



(2.14) 



(2.15) 



Here Sik is the Kronecker delta symbol. Within the el- 



liptic critical-state model (2.141, the magnitude of the 



critical current density J c draws an ellipse on the J±-J\\ 
plane. This model makes use of two phenomenological 
parameters, namely the extreme values J c ± and J c \\ 
for the radius of the ellipse drawn by the magnitude of 
the critical current density. In the numerical calcula- 



tions for solving the system of equations ( 2.7 1-( 2.10 1 



the relation (2.141 is rewritten in the form 



E i = E(J)(J- 1 ).,J k 



(2.16) 



E(J) = 



(2.17) 



0, J<J c (B,<j>) 
p(J - J c {B,<t>)), J>J c (B,<f>) 

where ( J~ j is the inverse of the matrix (J c )ik in 
(2.141. The magnitude of the critical current density, 
J c (B,(j>), is given by the expression 



J C (B,, 



cos 2 (</>) sin 2 (cf>) 



n-l/2 



(2.18) 



Here, <f> denotes the angle of the critical current den- 
sity J with respect to the direction of the flux den- 
sity B. If J c _l = J c |i, the elliptic critical-state model 



(2.161 goes over into the Bean's (circular) critical-state 
model (2.12). Besides, the calculations of electromag- 



netic fields with J close to J c are also independent of 
the auxiliary parameter p in Eq. (2.171. 



C. Extended elliptic model 

The elliptic critical-state model, described in previous 
subsection, has recently been extended in Refs. |28ll3"T] 
by introducing the general relations 



E± = 

En = 



P±J±, 
P\\ J \\> 



(2.19) 
(2.20) 



where p± and p\\ are nonlinear effective resistivities, 
having a ratio r = p\\/p± independent of J just above 
J c as it was experimentally found [55] . A model for the 
effective resistivities is given by [3T] 



Ej = 



Jii | -J cc )sign(J||), 



0, <| Jj_ |< J cd 

Pd(\ J± I -Jcd)sign(Ji), | J± |> J c d 

(2.21) 

0<| J|| |< Jcc 
I ^|| |> Jcc 

(2.22) 

Here, the subscripts "d" and "c" respectively refer to de- 
pinning and cutting. Besides, J c d = J c (B,(j)) | sin(^) | 
and J cc = J c {B,<f) \ cos(0) |, where J C (B,4>) is de- 
fined according to the elliptic critical-state model as in 
Eq. ( |2.18[ ). If | J - J c | /J c < 1, the extended ellip- 
tic critical-state model reduces to the original one [Eqs 
( 2.16| and ( 2.17 1] by replacing pa and p c in Eqs. (2.21 1 
and (2.22) with pJ c /J c ± and pJ c /J c \\, correspondingly. 
Hence, in the case of the original elliptic model, the 
ratio r = p\\/p± at J > J c is equal to J c ±/J c \\- On the 
other hand, the extended elliptic critical-state model 
is capable to modify the relation between the compo- 
nents of the electric field E and the current density J 
with the aid of the additional parameter r. 

D. Rectangular model 

The generalized double critical-state model (GDCSM) 
[H El 13 H] uses two phenomenological parameters, 
namely the critical values, J c \\ and J c j_, of the electri- 
cal current density along and perpendicular to the local 
magnetic induction. Within this model, each compo- 
nent of the electrical current density is determined by 
its own electric field as 



J± - 
J\\ 



J C _L sign (E^), 
= J c \\ sign(Si ). 



(2.23) 
(2.24) 



Evidently, the magnitude of the critical current density 
traces a rectangle in the J±_-J\\ plane. The parameter 
J c i_ determines the threshold for depinning of vortices, 
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whereas J c m indicates the onset of flux-line cutting in 
the vortex array. In calculating the electromagnetic 



fields within the GDCSM, the material equation (2.23) 
is written in the form 



E 



En = 



J\\ I -^ c ||)sign(J|| 



0, 0<|J ± |<J c± 
P±(\ J± I -J c x)sign{J±_), | J x \> Jc± 

(2.25) 

0<| J\\ \< J c \\ 

I J\\ \> Jc\\ 

(2.26) 

The quantities p± and pu are effective flux-flow and 
flux-line-cutting resistivities of the material. However, 
unlike the above-commented critical-state models, the 
GDCSM allows the existence of zones in the J±-J\\ 
plane where either flux cutting or flux transport ex- 
clusively occur. The latter is possible due to the as- 
sumption of the GDCSM that the threshold for flux 
depinning, J c j_ (flux cutting, J c m) is independent of the 



component Jy ( Jj_) [compare Eqs. (2.25) and (2.26) 
with Eqs. ( 2.21[ ) and (2.22) where J c d ana J cc depend 
on the angle <fi — arctan(Ji/Jii)]. 

III. Numerical results and comparison 

WITH EXPERIMENT 

In the present section we will apply the flux-line- 
cutting critical-state models, commented above, to ex- 
plain experimental magnetization curves |33| of a PbBi 
superconducting disk, rotating in the presence of an ex- 
ternal magnetic field H a , which is oriented parallel to 
the disk plane (along the z-axis) and perpendicular to 
the axis of rotation. 



v 
=1° 



Experiment 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 



^o H z CD 



Model 



0.06 



0.03 



0.00 



-0.03 



-0.06 




0.00 0.02 0.04 



0.06 0.08 



0.10 0.12 



0.14 



A. Experimental results 

Fig. [IJl exhibits a standard magnetization curve, which 
was measured in Ref. [33], for a PbBi disk of thick- 
ness d = 0.8mm. The hysteresis in Fig. [T]i clearly 
corresponds to the magnetization curve of a type-II ir- 
reversible superconductor since its return crosses over 
and remains in the paramagnetic region as a result of 
the strong flux pinning. In the experiment, the isotropy 
of the PbBi disk was also verified by comparing stan- 
dard magnetization curves with H a directed along dif- 
ferent diameters of the disk. 

Panels (a)-(c) in Fig. [2] show graphs of the mag- 
netization components, < M y >=< B y > //io and 
— < M z >= H a — < B z > //io, versus the angle 8 
of rotation, measured in the work [33 for the PbBi 
disk, rotating in the magnetic field H a . The measure- 
ments started in the nonmagnetic initial state which is 
reached after cooling the superconductor at the fields 
H a /Hp = 0.5 (panel a), 1.0 (panel b), and 2.0 (panel 
c), where Hp (poHp = 0.1015T |33|) is the penetra- 
tion field. The initial state is supposed to be non- 
magnetic because no Meissner effect (flux expulsion) 
was observed after field cooling, within the accuracy 
(A < M >< 1 Gauss) of the experiment. 



Figure 1: Standard magnetization curves (a) for a PbBi 
disk, taken from Ref. [33]. Theoretical magnetiza- 
tion curves (b) obtained with a critical current density 



J C ±(B) as in Eq. (3.31 1 



As it is seen in Fig. [2] for the smallest value of H a 
(= 0.5-ffp, panel a), both magnetization components 
have a nonmonotonic behavior as functions of 9. Such 
a behavior of magnetization has also been observed in 
Ref. [TTJ during the initial rotation of a Nb disk under- 
going slow oscillations in a parallel field. The depen- 
dence of the magnetization on 9 radically changes at 
larger values of H a . So (see Fig. [2jb), at H a = H P the 
functions < M y > (6) and — < M z > (9) initially grow 
with 9 and later (at 9 > 150°) they practically become 
constants with close values (< M y >« — < M z >). 
Also note that M y has a maximum at 9 ps 75°. For 
H a larger than the penetration field Hp (panel c), the 
function — < M z > (9) takes values smaller than those 
for < My > (9). Both of them are almost constant 
functions, except at small rotation angles because of 
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their fast initial growth. Thus, the maximum of < 
My > is shifted to a smaller value of 6 (« 40°). 

B. Theoretical predictions 

The models described in the previous section can be 
applied to explain the experimental results (Fig. [2| 
if we fix the sample and rotate the external magnetic 
field H a (2.1 1 by an angle a s — — 9 instead of fixing 



the magnetic field and rotating the superconducting 



sample. 

- <M Z 
tities: 



Then, the experimental values < M y > and 
> should respectively correspond to the quan- 



< M„ >= 



< M, >= H a 



1 

Hod 
1 

H d 



dxB'(x), 



dxB z (x), 



(3.27) 



(3.28) 



where 
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Figure 2: PbBi rotational curves measured in Ref. [33J 



AI, 



B' y = oc s X x ■ B = B(x) sin[a(x) - a s ], (3.29) 

B' z = d s B = B(x) cos[a(x) - a s }. (3.30) 

The calculations of magnetization components < 
> and — < M z > with the critical-state models, 
discussed in Sec. |IL] require the employment of the pa- 
rameters J C ±(B) and J c n(B), depending on the mag- 
netic induction. The former, J C ±(B), is determined 
from the experimental curves of magnetization versus 
the applied field, varying along one direction only as 
in Fig. [T] (In this case, flux cutting does not occur and, 
consequently, the depinning effects are completely re- 
sponsible for the magnetic response of the supercon- 
ductor). The standard magnetization curves are well 
reproduced by any one of the critical-state models (see 
above) with 



Jc±(B) = 



(1 + B/hoHpY 



(3.31) 



Jc_l(0) = 47.11 x 10 7 A/m 2 , and n ± = 2 (compare 
panels (a) and (b) of Fig. [T]). Other parameters of the 
critical state models are found by adjusting theoretical 
magnetization curves to the experimental ones (Fig. 

1. Circular model 



Within the Bean's circular critical-state model (2.11), 
there is only one phenomenological parameter, i.e. 
J C (B) = Jc±(B) = J C \\{B). Then, J C (B) has the form 



(3.31 1 with the same values for the parameters J c _l(0), 
and n±. 

Fig. [3] shows our numerical results for < M y > 
and — < M z >, obtained with the Bean critical-state 
model. At first glance, it seems that the circular model 
qualitatively reproduces the experimental magnetiza- 
tion curves (Fig. [2]). However, there are important 
differences between its predictions and the experiment. 
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Thus, for example, the "oscillations" of the magnetiza- 
tion components (Fig. [3ja) have small amplitudes com- 
pared with the experimental ones. Besides, at H a = 
O.SHp the functions < M y > (9) and - < M z > (9) 
approximate each other but at relatively large rotation 
angles 9 > 300° . Finally, when the applied field has an 
amplitude larger than H p (see panel c), the magneti- 
zation components are rather small in magnitude and 
their initial growth, before the saturation, occurs in a 
very small interval of 9 (< 20°). 



2. Elliptic model 

The calculations of magnetization components < M y > 
and — < M z > within the elliptic flux-line-cutting 
critical-state model ( |2.14 1 are shown in Fig. |4j Here, 
we used the same J C ±(B) as in Eq. (3.31 1 and J C \\(B) 
of the form 
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Figure 3: Curves of the average magnetization compo- 
nents versus the rotation angle, calculated with Bean's 
critical-state model. 



^||(0) 



(1 + B/ f i H P y 



(3.32) 



with J c ||(0) = 1.5J c j_(0) and nil = 1. This choice pro- 
vides a good agreement between experimental (Fig. |2| 
and theoretical (Fig. |i| curves. Thanks to the use of 
a second parameter (J c j|), the elliptic model is able to 
generate the "oscillations" of the magnetization com- 
ponents (Fig. |4]a) with amplitude close to that ob- 
served in the experiment (panel (a) in Fig. [2J). Notice 
that < My > and — < M z > approach each other 
at 9 > 150° with Ho = l.ObHp in good concordance 
with the measurements (see Fig. |2]b, corresponding 
to H = Hp). In addition, when H = 2.0H p (panel 
(c) in Fig. |2J, the difference between < M y > and 



- < M z > at 6 > 45° 
(Fig. gJ C ). 



is as large as in the experiment 



3. Extended elliptic model 

As was commented in Sec. |IL] both elliptic and cir- 
cular critical-state models are particular cases of the 
extended elliptic one. Therefore, the results presented 
in Fig. [3j predicted by the circular model, can also 
be calculated by using the new model (Eqs. (2.211 
and ( 2.22 1 ) with J c ± = J c ii as in E< 
r = p\\/p± = Pel Pd being equal to one (r 
J > J c . The condition r — 1 guarantees that the elec- 
tric field E and current density J be parallel as it is 



(3.311 and 
1) at 



postulated by Bean's critical-state model (2.11l. In 
addition, graphs in Fig. [4] (original elliptic model pre- 
dictions) , which quantitatively reproduce experimental 
measurements (Fig. [2]), are also obtained with the ex- 
tended elliptic critical-state model (Eqs. (2.21) and 
( |2.22| ) if r = J c ±/J c \\ (i-e. p c /pd = J c x/J c \\)- Accord- 
ing to the parameters J C ±(B) (3.311 and J C \\(B) (3.32), 
used for calculating magnetization curves in Fig. 4j the 
ratio r is here smaller than 1 (r < 1). 

It is interesting to study the effect of the parameter 
r, controlling the relation between the electric field E 



G 



Elliptic model 



Extended elliptic model (r=1) 




Elliptic model 



Extended elliptic model (r=1) 



o 



CD 

co 



0.8- 


b 




| i i | i i | i i 

H a /H p =1.05 - 


0.6- 


/<M> 

/ y 






0.4- 




-<M > 

Z 




0.2- 








0.0- 











45 90 135 180 225 270 315 360 

e (degrees) 



g 
'to 

N 

CD 
C 
D) 
CO 




45 90 135 180 225 270 315 360 

e (degrees) 



Elliptic model 



Extended elliptic model (r=1) 



0.8 



X, 0.6 

o 

00 0.4 
_N 

CD 

D) 0.2 
CO 



1 1 I 1 • 1 1 1 I 1 1 

c 


i ■ ■ i ■ ■ i i i i i i 
H a /H p =2.0 - 








' Q_ 

X 








■ / <M> 




o 


/ y 










_N 






' 

CD 


: // 




C 






D) 






CO 


\// 


-<M > 

z 




i i i i i i i i i i i i i i i i i i i i i i i 





45 90 135 180 225 270 315 360 

e (degrees) 



0.8 



0.0 -■ 



■ ■ 1 1 1 I 1 1 I 1 1 

c 


H a /H p =2.0 - 


/ <M > 

/ y 


'// 


- < M > 

Z 





45 90 135 180 225 270 315 360 

e (degrees) 



Figure 4: Curves of the average magnetization com- 
ponents versus the rotation angle, calculated with the 
original elliptic critical-state model. 



Figure 5: Curves of the average magnetization compo- 
nents versus the rotation angle, calculated with the ex- 
tended elliptic critical-state model using a ratio r = 1. 
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Figure 6: Curves of the average magnetization com- 
ponents versus the rotation angle, calculated with the 
generalized double critical-state model. 



and the current density J at J > J c . For this reason, 
we have calculated magnetization curves (Fig. |5| by 
applying the extended elliptic model with the same pa- 
rameters J C ±(B) and J C \\(B) as those employed in Fig. 
[4] but with the parameter r = p c /pd = 1. In other 
words, the magnetization curves in Fig. [5] correspond 
to an anisotropic critical-state model with J c ±/J c \\ < 1) 
but the parameter r = 1, indicating that E and J are 
parallel when J > J c . From the comparison of Fig. 
[5] with ^] we note that magnetization curves signifi- 
cantly depend upon the parameter r when the applied 
magnetic field is large enough (H a > Hp as in panels 

(b) and (c)). So, in order the magnetization compo- 
nents, < My > and — < M z >, to have the same 
value at large angles of rotation, the applied magnetic 
field H a for r = 1 (Fig. [5]b) should be larger than the 
field used in Fig. HI b. Besides, the value of < M y > 
and — < M z > (w OAHp), at sufficiently large an- 
gles 6, turns out to be smaller than that (« Q.bHp) 
predicted by the original elliptic model (Fig. |4]b) . 
At H a = 2.0-ffp, there is also a noticeable difference 
between magnetization ^-components (compare panels 

(c) of Figs. [f]and[5]). 

4- Rectangular model 

For completeness of our study, we have employed the 
GDCSM (rectangular model) , which also uses two crit- 
ical current densities, namely J C ±(B) and J C \\(B). The 
former is determined from the curves of magnetization 
versus the applied field, varying along one direction 
only (Fig. [TJ. In our case, the magnetic dependence of 
J c ± is the same as in Eq. (3.311. To reproduce the main 
features of the experiment (Fig. [2]), the other parame- 
ter is chosen as in Eq. ( pT32) , but Jcm(0) = 1.32J c _l(0) 
andrt|| = 1.06 (compare Figs. [2]and6j). Although these 
values are different from those usedwithin the elliptic 
critical-state model, the parallel critical current den- 
sity J c || remains being larger than the perpendicular 
one J c ±. It should be noted that the GDCSM predicts 
the equality of < M y > and - < M z > (w 0.5H P ) 
with an external field H a = 1.29-ffp > Hp at rela- 
tively large rotation angles 9 > 270° (see Fig. |6]b), 
in contrast to the experiment where such a behavior 
occurs from 9 s» 150°. Besides, the numerical calcu- 
lations for H a = 0.5Hp (panel (a) in Fig. [6]) had to 
be stopped at 9 ~ 338° because the solution further 
diverged. 

C. Magnetic induction profiles 

The fact that the elliptic critical-state model is able to 
quantitatively reproduce the experiment, with the use 
of a parallel critical current density J c y (B) larger than 
the perpendicular one J c ± (B) , illustrates how flux- line 
cutting influences on the magnetic behavior of a rotat- 
ing superconductor. To explain the features observed 
in both experimental (Fig. [2]) and theoretical (Fig. [1]) 
magnetization curves, we shall analyze the evolution 
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Figure 7: Profiles of the angle a (panel a), magnitude B (panel b) and components B' y (Eq. (3.29), panel c) 
and B' z (Eq. ( |3.30 1, panel d) of the magnetic induction, calculated with the original elliptic critical-state model 
at H a = O.bHp. 



of the profiles for the magnitude of the magnetic in- 
duction B(x), the tilt angle a{x), and the components 
B' y (x) ( |3.29| and g^(a:)( |3.30| , calculated within the 
original elliptic flux-line-cutting critical-state model 
(Figs. 09j. 



The calculated profiles of the magnetic induction 
in the case when the external magnetic field H a has 
a magnitude smaller than the penetration field H p 
(H a = 0.5H p ) are shown in Fig. [7j As the angle of ro- 
tation is increased, two [/-shaped minima in the B(x) 
profile (panel b) appear because of the flux consump- 
tion (decrement of B) which results from flux- line cut- 
ting [4] . The absolute value of the tilt angle a increases 
with 9 in the near-surface intervals < x < x m \ and 
x 7n 2 < x < d. However, in the intervals x m \ < x < x\ 
and X2 < x < x m 2, where there is flux consumption, 
the angle a is slightly modified. In the central inter- 
val, x\ < x < X2, neither B or a is altered. When 
9 w 360°, the minimum values of B inside the su- 
perconducting disk tend to zero and, as follows from 



Eq. (2.8 1, the magnitude of the derivative da/dx con- 
siderably increases at such points. Besides, at x = x m \ 



and x — x m 2 with B(x m \) = B{x m 2) ~ 0, the ac- 
curacy of our calculations is low and, therefore, the 
values —a(x m \) and —a{x m 2) turned out to be ap- 
parently higher than they should be (see curve 8 for 
9 = 360° in Fig. [7]a). The component B' z of the mag- 
netic induction, parallel to the applied magnetic field 
H a , decreases near sample surfaces because of the flux 
consumption (Fig. |7]d). Nevertheless, the most impor- 
tant change occurs in the central part of the sample 
(in Xi < x < X2) because of the sample rotation. So, 
at 9 = 180° (curve 4) the component B' z varies from 
B' z = /JoH a at the surfaces x — and x = d to the 
opposite value B' z = —[ioH a in the central region of 
the sample. When an entire cycle is finished, B' z again 
takes the value B' z = ^H a in the middle of the disk 
(curve 8). This cyclic behavior of B' z is responsible 
for the "oscillations" of the magnetization component 
< M z > (9) (panels (a) in Figs. [2] and HJ, being neg- 
ative for any value of the angle of rotation 9 > be- 
cause B' z < fioH a near surfaces, i.e. in the intervals 
< x < x\ and X2 < x < d. The component B' y also 
oscillates in the middle of the sample as 9 is increased 
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Figure 8: Profiles of the angle a (panel a), magnitude B (panel b) and components B' (Eq. (3.29), panel c) 



and B' z (Eq. ( |3.30 1, panel d) of the magnetic induction, calculated with the original elliptic critical-state model 
at H a = imHp. 



(Fig. [7]c). Such a behavior of B' y makes the magne- 
tization y-component < M y > oscillate with 9 (Figs. 
[2] a and |4] a). As it is seen in Fig. [7j;, there is an in- 
crement of B' y in the near-surface regions, producing a 
small positive value for < M y 
cycle, i.e. at 9 — 360° (see Figs. 



37271 after a complete 
^2ja and|4]a). 
Fig. [8] exhibits profiles calculated within the ellip- 
tic critical-state model for H a = 1.05H P . Due to the 



decrease of the critical current densities J c j_ (3.31 1 and 
J c ii (3.321 with the magnitude B of the magnetic in- 



duction, the slopes of the critical profiles for B(x) and 
a(x) near surfaces are smaller than the slopes observed 
in the corresponding profiles of Fig. [7] Therefore, the 
central region with unaltered B and a (see curves 1 
in panels (a) and (b) of Fig. |8| rapidly disappears as 
the rotation angle 9 is increased (see curves 2 therein). 
Also, the [/-shaped minima of B(x) coalesce forming a 
unique minimum at the center of the disk. The result- 
ing critical profile B(x) does not further change despite 
the fact that the disk continues rotating (see curves 5- 
8 in panel (b)). In this case, B' z (x) initially decreases 
(curves 1-2 in Fig. |8]d) inside the sample as 9 varies 
until it reaches the critical profile (curves 3-8). Hence, 
the dependence < M z > (9) has a monotonic behavior 



at 9 > 120° (see panels (b) in Figs. [2] and [i}. On the 
other hand, B' y (x) increases so that a huge maximum 
in the dependence < M y > (9) (Figs.[2]b and[I]b) ap- 
pears at 9 70°. At large rotation angles (9 > 180°), 
the profile B' y {x) becomes stationary and < M y > (9) 
is, practically, a constant function, having a value close 
to — < M z >. So, the magnitude of the magnetization, 
I < M > I , is independent of 9 when the rotation angle 
is sufficiently large. 

The profiles for the case when the external magnetic 
field is large enough, in comparison with the penetra- 
tion field Hp (as in Fig.|9|, have an evolution similar to 
that presented in Fig. [8J However, the central regions 
of unaltered magnetic induction rapidly disappear as 9 
is increased, in comparison with the results of Fig. [8] 
This fact is due to noticeable reduction of the critical 
current densities J c ± and J c n with B. 



IV. 



Conclusion 



We have applied the circular, elliptic, extended-elliptic, 
and rectangular critical-state models to study the mag- 
netic behavior of irreversible type-II superconductors 
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Figure 9: Profiles of the angle a (panel a), magnitude B (panel b) and components B' y (Eq. (3.29), panel c) 
and B' z (Eq. ( |3.30 1, panel d) of the magnetic induction, calculated with the original elliptic critical-state model 
at H a = 2.0H P . 



in a parallel rotating magnetic field. The numerical 
method employed here is based on the substitution of 
the vertical law, relating the electric field E and the 
current density J, for a nonlinear material equation 
having effective flux-cutting and flux-flow resistivities 
in the dissipative region. The substitution is justified 
when the applied magnetic field H a slowly varies either 
in magnitude or direction, inducing electric fields of 
sufficiently small magnitude inside the superconductor. 
Within the elliptic (circular) critical-state model such 
resistivities are not independent of each other and have 
a ratio r = p\\/p± equal to J c a_/J c \\ (=1 f° r the circu- 
lar model) at J just above its critical value J c . On the 
other hand, within the extended elliptic critical-state 
model the ratio r is an independent parameter to be 
determined. The rectangular critical-state model also 
uses two independent resistivities, pii and p±. However, 
unlike the other critical-state models, the GDCSM as- 
sumes that flux cutting and flux depinning do not affect 
each other. 

The comparison of the predictions of the mentioned 
critical-state models with experimental measurements 
of magnetization for a rotating PbBi disk in a fixed 
magnetic field [33 shows that the original critical-state 
model can reproduce the main features of the magne- 



tization curves. The circular and rectangular critical- 
state models only achieve a qualitative description of 
the experiment. The extended elliptic model, being 
more general than the original elliptic one, has allowed 
us to study the effect of the relation between E and J 
in the dissipative region. However, additional theoret- 
ical and experimental studies are needed to elucidate 
on the effects associated with both flux-cutting and 
flux-flow resistivities. 
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